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Abstract 

The use of Feshbach resonances for tuning the interparticle interaction 
in ultracold Fermi gases has led to remarkable developments, in particular 
to the creation and Bose-Einstein condensation of weakly bound diatomic 
molecules of fermionic atoms. These are the largest diatomic molecules ob- 
tained so far, with a size of the order of thousands of angstroms. They rep- 
resent novel composite bosons, which exhibit features of Fermi statistics 
at short intermolecular distances. Being highly excited, these molecules 
are remarkably stable with respect to coUisional relaxation, which is a 
consequence of the Pauli exclusion principle for identical fermionic atoms. 
The purpose of this review is to introduce theoretical approaches and de- 
scribe the physics of molecular regimes in two-component Fermi gases and 
Fermi-Fermi mixtures, focusing attention on quantum statistical effects. 
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1 Introduction 

1.1 State of the art 

The field of quantum gases is rapidly expanding in the direction of ultracold 
clouds of fermionic atoms, with the goal of revealing novel macroscopic quan- 
tum states and achieving various regimes of superfluidity. The initial idea was 
to achieve the Bardeen-Cooper-Schrieffer (BCS) superfluid phase transition in 
a two-component Fermi gas, which requires attractive interactions between the 
atoms of different components. Then, in the simplest version of this transition, 
at sufficiently low temperatures fermions belonging to different components and 
with opposite momenta on the Fermi surface form correlated (Cooper) pairs 
in the momentum space. This leads to the appearance of a gap in the single- 
particle excitation spectrum and to the phenomenon of superfluidity (see, for 
example, [T]). In a dilute ultracold two-component Fermi gas, most efficient is 
the formation of Cooper pairs due to the attractive intercomponent interaction 
in the s-wave channel (negative s-wave scattering length a). However, for typi- 
cal values of a, the superfluid transition temperature is extremely low. For this 
reason, the efforts of many experimental groups have been focused on modify- 
ing the intercomponent interaction using Feshbach resonances. The scattering 
length a near a Feshbach resonance can be tuned from — oo to -|-oo. This has 
led to exciting developments (see [2] for review), such as the direct observation 
of superfluid behavior in the strongly interacting regime (n|ap > 1, where n is 
the gas density) through vortex formation [3], and the study of the influence 
of imbalance between the two components of the Fermi gas on superfluidity 

We focus here on the remarkable physics of weakly bound diatomic molecules 
of fermionic atoms. This initially unexpected physics connects molecular and 
condensed matter physics. The weakly bound molecules are formed on the posi- 
tive side of the resonance (a > 0) O [TOl [HI [12] and they are the largest diatomic 
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molecules obtained so far. Their size is of the order of a and it reaches thou- 
sands of angstroms in current experiments. Accordingly, their binding energy 
is exceedingly small (10 /iK or less). Being composite bosons, these molecules 
obey Bose statistics, and they have been Bose-condensed in experiments with 
4°K2 [Hm] at JILA and with ^Liz at Innsbruck [THIIlj, MIT [HlHj, ENS 
[T^ . Rice HO], and Duke [21]. Nevertheless, some of the interaction properties 
of these molecules reflect Fermi statistics of the individual atoms forming the 
molecule. In particular, these molecules are found to be remarkably stable with 
respect to collisional decay. Being in the highest rovibrational state, they do not 
undergo collisional relaxation to deeply bound states on a time scale exceeding 
seconds at densities of about 10^'^ cm~'^. This is more than four orders of mag- 
nitude longer than the life time of similar molecules consisting of bosonic atoms. 
The key idea of our discussion of homonuclear diatomic molecules formed in a 
two-component Fermi gas by atoms in different internal (hyperfine) states is to 
show how one obtains an exact universal result for the elastic interaction be- 
tween such weakly bound molecules and how the Fermi statistics for the atoms 
provides a strong suppression of their collisional relaxation into deep bound 
states. It is emphasized that the repulsive character of the elastic intermolecu- 
lar interaction and remarkable collisional stability of the molecules are the main 
factors allowing for their Bose-Einstein condensation and for prospects related 
to interesting manipulations with these molecular condensates. 

Currently, a new generation of experiments is being developed for studying 
degenerate mixtures of different fermionic atoms [22l [23] , with the idea of re- 
vealing the influence of the mass difference on superfluid properties and finding 
novel types of superfluid pairing. On the positive side of the resonance one 
expects the formation of heteronuclear weakly bound molecules, which attracts 
a great deal of interest, in particular for creating dipolar gases. We present 
an analysis of how the mass ratio for constituent atoms influences the elastic 
interaction between the molecules and their collisional stability. The discussion 
is focused on molecules of heavy and light fermions, where one expects the for- 
mation of trimer bound states and the manifestation of the Efimov effect. We 
then show that a many-body system of such molecules can exhibit a gas-crystal 
quantum transition. Remarkably, the atomic system itself remains dilute, and 
the crystalline ordering is due to a relatively long-range interaction between the 
molecules originating from exchange of light fermions. Realization of the crys- 
talline phase requires a very large mass ratio for the atoms forming a molecule 
in order to suppress the molecular kinetic energy. This can be achieved in an 
optical lattice for heavy atoms, where the crystalline phase of a dilute molecular 
system emerges as a superlattice, and we discuss the related physics. 

The chapter is concluded by an overview of prospects for manipulations with 
the weakly bound molecules of fermionic atoms. The leading ideas include the 
achievement of ultra-low temperatures and BCS transition for atomic fermions, 
creation of dipolar quantum gases, as well as observation of peculiar trimer 
bound states in an optical lattice. 
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Figure 1: The dependence of the scattering length on the magnetic field near 
a Feshbach resonance. The symbols I, II, and III label the regime of a weakly 
interacting degenerate atomic Fermi gas, strongly interacting regime of BCS- 
BEC crossover, and the regime of weakly bound molecules. At sufficiently low 
temperatures region I corresponds to the BCS superfluid pairing, and region III 
to Bose-Einstein condensation of molecules. 

1.2 Feshbach resonances and diatomic molecules 

At ultralow temperatures, when the de Broglie wavelength of atoms greatly ex- 
ceeds the characteristic radius of interatomic interaction forces, atomic collisions 
and interactions are generally determined by the s-wave scattering. Therefore, 
in two-component Fermi gases one may consider only the interaction between 
atoms of different components, which can be tuned by using Feshbach reso- 
nances. 

The description of a many-body system near a Feshbach resonance requires 
a detailed knowledge of the 2-body problem. In the vicinity of the resonance, 
the energy of a colliding pair of atoms in the open channel is close to the 
energy of a molecular state in another hyperfine domain (closed channel). The 
coupling between these channels leads to a resonant dependence of the scattering 
amplitude on the detuning 5 of the closed channel state from the threshold of 
the open channel, which can be controlled by an external magnetic (or laser) 
field. Thus, the scattering length becomes field dependent (see Fig. 1). 

The Feshbach effect is a two-channel problem which can be described in 
terms of the Breit-Wigner scattering [2ll [25] , and various aspects of such prob- 
lems have been discussed by Feshbach [26] and Fano [27]. In cold atom physics 
the idea of Feshbach resonances was introduced in Ref. [28] , and optically in- 
duced resonances have been analyzed in Refs. [29l [30l [Ml [32] . 

At resonance the scattering length changes from -|-c» to —oo, and in the 
vicinity of the resonance one has the inequality n|ap > 1 where n is the gas 
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density. The gas is ten said to be in the strongly interacting regime. It is still 
dilute in the sense that the mean interparticle separation greatly exceeds the 
characteristic radius of the interparticle interaction R^. However, the amplitude 
of binary interactions (scattering length) is larger than the mean separation 
between particles, and in the quantum degenerate regime the conventional mean 
field approach is no longer valid. 

For large detuning from resonance the gas is in the weakly interacting regime, 
i.e. the inequality n\a\^ ^ 1 is satisfied. On the negative side of the resonance 
(a < 0), at sufficiently low temperatures of the two-species Fermi gas one expects 
the BCS pairing between distinguishable fermions, well described in literature 
[T] . On the positive side (a > 0) two fermions belonging to different components 
form diatomic molecules. For a ^ Re these molecules are weakly bound and 
their size is of the order of a. 

The crossover from the BCS to EEC behavior has recently attracted now a 
great deal of interest, in particular with respect to the nature of superfluid pair- 
ing, transition temperature, and elementary excitations. This type of crossover 
has been earlier discussed in the literature in the context of superconductivity 
[33l [34l [35l [36] and in relation to superfluidity in two-dimensional films of ^He 
[571 [55] . The idea of resonant coupling through a Feshbach resonance for achiev- 
ing a superfluid phase transition in ultracold two-component Fermi gases has 
been proposed in Refs. [351110]) and for the two-dimensional case it has been 
discussed in Ref. [H] . 

The two-body physics of the Feshbach resonance is the most transparent if 
the (small) background scattering length is neglected. Then for low collision 
energies e the scattering amplitude is given by [25] : 

Fis)=- "y^^ , (1) 

where the quantity h'^ / ^/Ijl = W characterizes the coupling between the open 
and closed channels and /x is the reduced mass of the two atoms. The scatter- 
ing length is a = — F(0). In Eq. ([T|) the detuning 6 is positive if the bound 
molecular state is below the continuum of the colliding atoms. Then for (5 > 
the scattering length is positive, and for (5 < it is negative. Introducing a 
characteristic length 

R* = h^/2fiW (2) 

and expressing the scattering amplitude through the relative momentum of par- 
ticles k = ^J2^£/h, we can rewrite Eq. ([T]) in the form: 

The validity of Eq. does not require the condition kR* <C 1. At the same 
time, this equation formally coincides with the amplitude of scattering of slow 
particles by a potential with the same scattering length a and an effective range 
R = —2R*, obtained under the condition k\R\ ^ 1. 
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The length R* is an intrinsic parameter of a Feshbach resonance. It charac- 
terizes the width of the resonance. From Eqs. H]) and ([2]) we see that small W 
and, consequently, large R* correspond to narrow resonances, whereas large W 
and small R* lead to wide resonances. The term "wide" is generally used when 
the length R* drops out of the problem, which according to Eq. ^ requires 
the condition kR* <C 1. In a quantum degenerate atomic Fermi gas the charac- 
teristic momentum of particles is the Fermi momentum kp = (Stt^tt,)^/^. Thus, 
in the strongly interacting regime and on the negative side of the resonance 
(a < 0) , for a given R* the condition of the wide resonance depends on the gas 
density n and takes the form kpR* < 1 [i^ H51 HI H5 1 HS] . 

For a > one has weakly bound molecular states (it is certainly assumed 
that the characteristic radius of interaction R^ <C a), and for such molecular 
systems the criterion of the wide resonance is different [47l |49] . The binding 
energy of the weakly bound molecule state is determined by the pole of the 
scattering amplitude ^ . One then finds [47l |49] that this state exists only for 
a > 0, and under the condition 

i?* < a (4) 

the binding energy is given by 

£0 = h^/2^ia^. (5) 

The wavefunction of such weakly bound molecular state has only a small admix- 
ture of the closed channel, and the size of the molecule is ~ a. The characteristic 
momenta of the atoms in the molecule are of the order of a^^ and in this respect 
the inequality ([4]) represents the criterion of a wide resonance for the molecular 
system. 

Under these conditions atom-molecule and molecule-molecule interactions 
are determined by a single parameter - the atom-atom scattering length a. In 
this sense, the problem becomes universal. It is equivalent to the interaction 
problem for the two-body potential which is characterized by a large positive 
scattering length a and has a potential well with a weakly bound molecular state. 
The picture remains the same when the background scattering length can not be 
neglected, although the condition of a wide resonance can be somewhat modified 

m- 

Most ongoing experiments with Fermi gases of atoms in two different inter- 
nal (hyperfine) states use wide Feshbach resonances [51]. For example, weakly 
bound molecules ^Li2 and have been produced in experiments [101 El 

[l2l[13l[14l[17l[T8l[l5l[l6l[19l[2n| by using Feshbach resonances with a length 
R* < 20A, and for the achieved values of the scattering length a (from 500 to 
2OOOA) the ratio R* /a was smaller than 0.1. In this review we will consider the 
case of a wide Feshbach resonance. 
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2 Homonuclear diatomic molecules in Fermi gases 



2.1 Weakly interacting gas of bosonic molecules. Molecule- 
molecule elastic interaction 

As we have shown in the previous section, the size of weakly bound bosonic 
molecules formed at a positive atom-atom scattering length a in a two-species 
Fermi gas (region III in Fig. 1) is of the order of a. Therefore, at densities such 
that na^ ^ 1, the atoms form a weakly interacting gas of these molecules. More- 
over, under this condition at temperatures sufficiently lower than the molecu- 
lar binding energy Eq and for equal concentrations of the two atomic compo- 
nents, practically all atoms are converted into molecules [52]. This is defi- 
nitely the case at temperatures below the temperature of quantum degeneracy 
Td = 2nh?n'^/^ /M (the lowest one in the case of fermionic atoms with different 
masses, with M being the mass of the heaviest atom). One can clearly see this 
by comparing Td with Eq given by Eq ([5]). Thus, one has a weakly interacting 
molecular Bose gas and the first question is related to the elastic interaction 
between the molecules. 

For a weakly interacting gas the interaction energy in the system is equal 
to the sum of pair interactions and the energy per particle is ng {2ng for a 
non-condensed Bose gas), with g being the coupling constant. In our case this 
coupling constant is given hy g = ATrh?add/ {M + 171), where add is the scattering 
length for the molecule-molecule (dimer-dimer) elastic s-wave scattering, and 
M , m are the masses of heavy and light atoms, respectively. The value of add 
is important for evaporative cooling of the molecular gas to the regime of Bose- 
Einstein condensation and for the stability of the condensate. The Bose-Einstein 
condensate is stable for repulsive intermolecular interaction {add > 0), and for 
add < it collapses. 

We thus see that for analyzing macroscopic properties of the molecular Bose 
gas one should first solve the problem of elastic interaction (scattering) between 
two molecules. In this section we present the exact solution of this problem 
for homonuclear molecules formed by fermionic atoms of different components 
(different internal states) in a two-component Fermi gas. The case of M 7^ 
m will be discussed in Section 3. The solution for M = m was obtained in 
Refs. [IHl [49] assuming that the atom-atom scattering length a greatly exceeds 
the characteristic radius of interatomic potential: 

a > Re. (6) 

Then, as in the case of the 3-body problem with fermions [551 [Ml [53] [53], the 
amplitude of elastic interaction is determined only by a and can be found in the 
zero-range approximation for the interatomic potential. 

This approach was introduced in the two-body physics by Bethe and Peierls 
[57] . The leading idea is to solve the equation for the free relative motion of two 
particles placing a boundary condition on the wavefunction ^ at a vanishing 
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Figure 2: Set of coordinates for the four-body problem. 



interparticle distance r: 



rip 



1 

a 



which can also be rewritten as 



ip (X (1/r — 1/a), r 



0. 



(7) 



(8) 



One then gets the correct expression for the wavefunction at distances r ^ Re. 
when a ^ Re, Eq. ([5]) correctly describes the wavefunction of weakly bound 
and continuum states even at distances much smaller than a. 

We now use the Bethe-Peierls approach for the problem of elastic molecule- 
molecule (dimer-dimer) scattering which is a 4-body problem described by the 
Schrodinger equation 



-V 



C/(ri) + C/(r2)+X]t/[(ri+r2±%/2R)/2]-£;}* = 0, (9) 
± 



where m is the atom mass. Labeling fermionic atoms in different internal states 
by the symbols t and |, the distance between two given t and | fermions is ri, 
and r2 is the distance between the other two. The distance between the centers 
of mass of these pairs is R/-\/2, and (ri -I- r2 ± ■\/2R)/2 are the separations 
between t and | fermions in the other two possible ti pairs (see Fig. 2). The 
total energy is E = —2eo + e, with e being the colhsion energy, and Eq = 
—h'^/ma'^ the binding energy of a dimer. The wavefunction ^ is symmetric 
with respect to the permutation of bosonic ti pairs and antisymmetric with 
respect to permutations of identical fermions: 



-r2T\/2R ri-r2 
2 ' ^/2 



(10) 
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For the weak binding of atoms in the molecule assuming that the 2-body 
scattering length satisfies the inequality ([6]), at all interatomic distances (even 
much smaller than a) except for very short separations of the order of or smaller 
than Re, the motion of atoms in the 4-body system is described by the free- 
particle Schrodinger equation 



„o „9 „o mE 



* = 0. (11) 



The correct description of this motion requires the 4-body wavefunction ^' to 
satisfy the Bethe-Peierls boundary condition for the vanishing distance in any 
pair of t and | fermions, i.e. for ri ^ 0, r2 ^ 0, and ri -f r2 ± V^R 0. Due 
to the symmetry condition (jlOp it is necessary to require a proper behavior of 
\[' only at one of these boundaries. For ri the boundary condition reads: 

*(ri,r2,R) ^ /(r2,R)(l/47rri - l/4^a). (12) 

The function /(r2, R) contains the information about the second pair of particles 
when the first two are on top of each other. 
In the ultracold limit, where 

fca<l, (13) 

the molecule-molecule scattering is dominated by the contribution of the s-wave 
channel. The inequality ([13]) is equivalent to e <C eo and, hence, the s-wave 
scattering can be analyzed from the solution of Eq. (fTTjl with E = —2eo < 0. 
For large R the corresponding wavefunction takes the form 

^ ^Mn)M'r2){l~V2add/R); i?>a, (14) 
where the wavefunction of a weakly bound molecule is given by 

Mr) = ^^^exp{~r/a). (15) 



Combining Eqs. p2)) and (|T4l) we obtain the asymptotic expression for / at 
large distances R: 

/(r2, R) » (2/r2a) exp (-r2/a)(l - V2add/R); i? > a. (16) 

In the case of s-wave scattering the function / depends only on three vari- 
ables: the absolute values of r2 and R, and the angle between them. We now 
derive and solve the equation for /. The value of the molecule-molecule scat- 
tering length add is then deduced from the behavior of / at large R determined 
by Eq. dUD. 

We first establish a general form of the wavefunction '5 satisfying Eq. pT|) . 
with the boundary condition (jl2p and symmetry relations pop . In our case the 
total energy E = —2fi^/ma^ < 0, and the Green function of Eq. pT|) reads 

G{X) = {2TT)-'>/\Xa/V2y^/^Ky2{V2X/a), (17) 
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where X ^ \S — S'\, /2{V2X / a) is the deacaying Bessel function, and S = 

{ri, r2, R} is a 9-component vector. Accordingly, 15*— S"| = (ri — r'l)^ + (r2 — r'2)^ + (R — R')^. 

The 4-body wavefunction 5* is regular everywhere except for vanishing distances 

between t and j fermions. Therefore, it can be expressed through G(|5 — S'\) 

with coordinates S' corresponding to a vanishing distance between "f and | 

fermions, i.e. for r'l 0, r'2 0, and (r'l + r'2 ± \/2R')/2 ^ 0. Thus, for the 

wavefunction 5* satisfying the symmetry relations (jlOp we have 



-G(|5 - S'+l) - GdS" - S'_|)j /i(r', R'), (18) 

where Si = {0,r',R'}, S2 = {r',0,-R'}, and S± = {r72 ± R7V2, r72 T 
R'/\/2 =F r'\/2}. The function is a properly symmetrized finite solution of 
Eq. pT|) . regular at any distances between the atoms. For E < 0, non-trivial 
solutions of this type do not exist and we have to put = 0. The function 
/i(r2,R) has to be determined by comparing 'i' in Eq. (fT8|) at ri 0, with the 
boundary condition (|12|) . 

Considering the limit ri we extract the leading terms on the right hand 
side of Eq. . These are the terms that behave as 1 /ri or remain finite in this 
limit. The last three terms in the square brackets in Eq. (jl8p provide a finite 
contribution 

j cPr'cPR' /i(r',R') [g(|52 - ^21) - G{\S2 - S+\) - G{\S2 - ^-|)] , (19) 

where S2 = {0, r2,R}. To find the contribution of the first term in the square 
brackets, we subtract and add an auxiliary quantity 



/i(r2,R) / G(|5-5i|)dV(i3i?' = M£ai^exp(-V2ri/a). 
J inri 



(20) 



The result of the subtraction yields a finite contribution which for ri — s- can 
be written as 

J d^r'd^R'[h{r',R') - h{r2,'R)]G{\S - Si\) 

J dVd3i?'[/j(r',R') - hir2,R)]Gi\S2 - ^i|); n ^ 0, (21) 

with the symbol P denoting the principal value of the integral over dr' (or dR'). 
A detailed derivation of Eq. pT|) and the proof that the integral in the second 
line of this equation is convergent are given in Ref. [49] . 

In the limit ri 0, the right hand side of Eq. ([20]) is equal to 



/i(r2, R)(l/47rri - \/2/47ra). (22) 
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We thus find that for ri ^ the wavefunction of Eq. ^TE^ takes the form 

M^(ri,r2,R)^ ^^'"^'^^ +TZ- n 0, (23) 

47rri 

where TZ is the sum of regular n-independent terms given by Eqs. (fT9|) and ((2T|) . 
and by the second term on the right hand side of Eq. ([^^ . Equation ([221) must 
coincide with Eq. p2p , and comparing the singular terms of these equations we 
find /i(r2,R) = /(r2,R). As the quantity TZ must coincide with the regular 
term of Eq. p^ . equal to — /(r2R)/47ra, we obtain the following equation for 
the function /: 

J dVd3i?'|G(|5-5i|)[/(r',R')-/(r,R)] + [G(|5-52|) 

- J2 G{\S - 5±|)] /(r', R')} = (V2 - l)/(r, R)/4™. (24) 

Here S = {0, r, R}, and we omitted the symbol of the principal value for the 
integral in the first line of Eq. . 

As we have already mentioned above, for s-wave scattering the function 
/(r, R) depends only on the absolute values of r and R and on the angle be- 
tween them. Thus, Eq. ([24]) is an integral equation for the function of three 
variables. In order to find the molecule-molecule scattering length, it is more 
convenient to transform Eq. (|24p into an equation for the momentum-space 
function /(k, p) = / d'^r(i'^i?/(r, R) exp(ik • r/a -I- «p • R/\/2a), which yields 
the following expression: 



/(k±(p'-p)/2,p')dV 



2 +p'2/2 +(k ± (p'- p)/2)2 + (k ± (p'+ p)/2)2 
/(k', -p) d3fc' 2^2(1 ^^2^ pV2)/(k, p) 



2 + fc'2 + A:2+p2/2 ^2 + fc2+p2/2 + i 



(25) 



By making the substitution /(k, p) = ((5(p) + g(k,p)/p'^)/{l + k^) we reduce 
Eq. ([25]) to an inhomogcneous equation for the function ^(k, p): 

1 27r2(l + fc2+p2/2)g(k,p) 



(l + fc2+p2/4)2_(kp)2 p2(l + fc2)(^2 + fc2+p2/2 + l) 

^^^f g(k±(p'-p)/2,p0dV 

2^J p'2(2+p/2/2 +(k±(p'-p)/2)2+(k±(p'+p)/2)2)(H-(k± (p'-p)/2)2) 

g(k',-p)d3fc' 



p2(2 + fc'2 + fc2 -|-p2/2)(l + fc'2) • 



(26) 



In the case of s-wave scattering the function g(k, p) depends on the absolute 
values of k and p and on the angle between these vectors. For p — + this func- 
tion tends to a finite value independent of k. As one can easily establish on the 
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basis of Eq. ([TB| and the definition of ^(k, p), the molecule-molecule scattering 
length is given by add — — 27r^a limp^o p). Numerical calculations from 
Eq. (pS)) give with 2% accuracy [S5] : 

add = 0.6a > 0. (27) 

This result was first obtained in Refs. [48l|49] from the direct numerical solution 
of Eq. ([M)) by fitting the obtained /(r, R) with the asymptotic form (|16p at large 
R. The calculations show the absence of 4-body weakly bound states, and the 
behavior of / at small R suggests a soft-core repulsion between dimers, with a 
range ~ a. 

The result of Eq. (P7)l is exact, and it indicates the stability of molecular BEC 
with respect to collapse. Compared to earlier studies which assumed add = 2a 
[36l [58] , Eq. (|27|) gives almost twice as small a sound velocity of the molecular 
condensate and a rate of elastic collisions smaller by an order of magnitude. 
The result of Eq. ([27]) has been confirmed by Monte Carlo calculations [59] and 
by calculations within the diagrammatic approach |60|, I61j . An approximate 
diagrammatic approach leading to a = 0.75a has been developed in Ref. [58j . 

2.2 Suppression of collisional relaxation 

Weakly bound dimers that we are considering are diatomic molecules in the high- 
est ro- vibrational state (see Fig. 3). They can undergo relaxation into deeply 
bound states in their collisions with each other: for example, one of the colliding 
molecules may relax to a deeply bound state while the other one dissociates |62j . 
The released energy is the binding energy of the final deep state, which is of the 
order of fi-^/mi?^. It is transformed into the kinetic energy of the particles in the 
outgoing collision channel and they escape from the trapped sample. Therefore, 
the process of collisional relaxation of weakly bound molecules determines the 
lifetime of a gas of these molecules and possibilities to Bose-condense this gas. 

We now show that collisional relaxation is suppressed due to Fermi statistics 
for atoms in combination with a large size of weakly bound molecules [IHl SH] ■ 
The binding energy of the molecules is Eq = /raa^ and their size is ~ a ^ i?e- 
The size of deeply bound states is of the order of Re- Therefore, the relaxation 
process may occur when at least three fermionic atoms are at distances ~ Re 
with respect to each other. As two of them are necessarily identical, due to 
the Pauli exclusion principle the relaxation probability acquires a small factor 
proportional to a power of {qRe), where q ^ 1/a is a characteristic momentum 
of the atoms in the weakly bound molecular state. 

Relying on the inequality a ^ R^ we outline a method that allows one 
to establish the dependence of the relaxation rate on the scattering length a, 
without going into a detailed analysis of the short-range behavior of the system. 
It is assumed that the amplitude of the inelastic relaxation process is much 
smaller than the amplitude of elastic scattering. Then the dependence of the 
relaxation rate on a is related only to the a-dependence of the initial-state 4- 
body wavefunction ^. We again consider the ultracold limit described by the 
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Figure 3: Interaction potential J7 as a function of the distance R between two 
distinguishable fermionic atoms. The dashed line shows the energy level of the 
weakly bound molecule, and the solid line the energy level of a deeply bound 
state. 




Figure 4: Configuration space contributing to the relaxation probability. 

condition (I13p . where the relaxation process is dominated by the contribution 
of the s-wave molecule-molecule scattering. 

The key point is that the relaxation process requires only three atoms to 
approach each other to short distances of the order of R^. The fourth particle 
can be far away from these three and, in this respect, does not participate in 
the relaxation process. This distance is of the order of the size of a molecule, 
which is ^ a ^ Re- We thus see that the configuration space contributing to 
the relaxation probability can be viewed as a system of three atoms at short 
distances Re from each other and a fourth atom separated from this system 
by a large distance ~ a (see Fig. 4). In this case the 4-body wavefunction 
decomposes into a product: 

* = 77(z)*(3)(p^f])^ (28) 

where ^'^■^^ is the wavefunction of the 3-fermion system, p and ft are the hyper- 
radius and the set of hyperangles for these fermions, z is the distance between 
their center of mass and the fourth atom. The wavefunction ri{z) describes the 
motion of this atom and is normalized to unity. Note that Eq. (|28p remains 
valid for any hyperradius p ^ \z\ ^ a. 

The transition to a deeply bound 2-body state occurs in the system of three 
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atoms and does not change the wavefunction of the fourth atom, r]{z). Therefore, 
averaging the transition probabihty over the motion of the fourth particle, the 
rate constant for inelastic relaxation in dimer-dimer collisions can be written as 

are; =a(3) J \j^{z)\^ z = a^^\ (29) 

where a*-'^^ is the rate constant for relaxation in the 3-atom system. 

At interatomic distances ~ Re where the relaxation occurs, as well as at all 
distances where the hyperradius p <^ a, the wavefunction ^'^'^^ is determined by 
the Schrodinger equation with zero energy and, hence, depends on the scattering 
length a only through a normalization coefficient : 

*(3)=A(a)V'; P<a, (30) 

where the function ip is independent of a. The probability of relaxation and, 
hence, the relaxation rate constant are proportional to l^fp at distances ~ R^. 
We thus have 

a,ei =a^'^ oc |A(a)p. (31) 

The goal then is to find the coefficient A{a), which determines the dependence 
of the relaxation rate on a. 

For this purpose it is sufficient to consider distances where R^ and 

Eq. pO|) is still valid. Then, using the zero-range approximation we find the 
coordinate dependence of the three-body wavefunction ^f'^^. The derivation is 
presented in Ref. J49] and the result is: 

p<a, (32) 

where $i/(f2) is a normalized function of hyperangles, and the coefficient v de- 
pends on the symmetry of ^f^^^. The a-dependence of the prefactor A{a) can be 
determined from the following scaling arguments for the 4-body problem. The 
scattering length a is the only length scale in our problem and we can measure 
all distances in units of a. Using two rescaled coordinates, p = ap' and z = az', 
we see that in Eqs. (|32|) and ((28)) becomes a function of p/a, multiplied 
by A{a)a'^~^. The wavefunction 7y(z) is normalized to unity and hence it is a 
function of z/a, multiplied by a~'^/^. Accordingly, the 4-body wavefunction ^ 
of Eq. (|28p is a function of rescaled coordinates, multiplied by the coefficient 
A{a)a''^^^^ . By applying the same rescaling to Eq. and using Eq. (fTCj) 
we see that the same coefficient should be proportional to a~^. Therefore, 
A{a) (X a"""^/^ and arei oc a~^, where s = 2;^ + 1. 

The strongest relaxation channel corresponds to the lowest value of ly. It 
is achieved in the case of p-wave symmetry in the three-body system described 
by the wavefunction ^'(^^ and is equal to v = 0.773, which leads to s = 2.55. 
Assuming that the short-range physics is characterized by the length scale Re 
and the energy scale /mR^ we can restore the dimensions and write: 

arei = C{hRe/m){Re/a)'; s = 2.55, (33) 
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where the coefRcient C depends on a particular system and can not be obtained 
using the zero-range approximation. 

Note that the p-wave symmetry in the three-body system corresponds to 
p-wave scattering of a fermionic atom of one of the molecules (referred to as 
the 3-rd fermion) on the other molecule. Then the 4-th particle also undergoes 
j)-wave scattering on this molecule in such a way that the total orbital angular 
momentum of the molecule-molecule collision is equal to zero. Since the 3-rd 
and 4-th fermions are bound to each other in the molecular state with a size 
~ a, the relative momentum of their collisions with the other molecule is ~ 1/a 
and such p-wave collisions are not at all suppressed. The relaxation channel 
corresponding to s-wave scattering of the 3-rd fermion on the molecule leads 
to = 1.1662 and hence to the relaxation rate proportional to a^^-^^ as in 
the case of ultracold atom-molecule collisions [48l|49]. Thus, for large a this 
mechanism can be omitted. The channels, where the 3-rd fermion (and the 4-th 
one) scatters on the molecule with orbital angular momentum / > 1, lead to 
even stronger decrease of the relaxation rate with increasing o and hence can 
be neglected. 

2.3 Collisional stability and molecular BEC 

Equation (|33p implies a remarkable collisional stability of weakly bound molecules 
consisting of fermionic atoms in two different internal states and a counterin- 
tuitive decrease of the relaxation rate with increasing the two-body scattering 
length a. For currently achieved values of the scattering length a ^ 1000 A, 
the suppression factor {Re /ay for the relaxation process is about 4 orders of 
magnitude. This effect is due to Fermi statistics for the atoms. It is not present 
for weakly bound molecules of bosonic atoms, even if they have the same large 
size. Indeed, as the size of weakly bound molecules is ~ a, identical fermionic 
atoms participating in the relaxation process have very small relative momenta 
k ^ 1/a. Hence, the probability that they approach each other to short dis- 
tances ~ Re where the relaxation transitions occur, should be suppressed as 
{kRe)^ ~ {Re/a^ compared to the case of molecules of bosonic atoms. The 
exponent s in Eq. ([55)) is different from 2 due to the Frank-Condon factor for 
the relaxation transition and three-body dynamics. 

The remarkable collisional stability of weakly bound molecules K2 and Li2 
consisting of two fermionic atoms has been observed in experiments at JILA 
[HIISlIIl], Innsbruck [TH [15l [16] , MIT [HI [18], ENS [TQl!T9l, Rice [20], and 
Duke [21]. At molecular densities n ~ 10"'^^ cm~^ the lifetime of the gas ranges 
from tens of milliseconds to tens of seconds, depending on the value of the 
scattering length a. A strong decrease of the relaxation rate with increasing a, 
following from Eq. p3p . is consistent with the experimental data. The potassium 
experiment at JILA [12] and lithium experiment at ENS [19] give the relaxation 
rate constant arei oc a~*, with s ~ 2.3 for K2, and s ~ 1.9 for Li2, in agreement 
with theory (s « 2.55) within experimental uncertainty. The experimental and 
theoretical results for potassium and lithium are shown in Fig. 5 and Fig 6. The 
absolute value of the rate constant for a ^Li2 condensate is arei ~ 1 x 10"-'^^ 
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Figure 5: Two-body decay rate of a '*°K2 ultracold molecular g function 
of the magnetic field detuning from the '^''K Feshbach resonance at 202 G. The 
dots indicate experimental values, and the solid line shows the theoretical results 
normalized to the experimental value at AB = —1.6 G. 
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Figure 6: Two-body decay rate a^ei of a ^Li2 molecular condensate as a function 
of interatomic scattering length near the ^Li Feshbach resonance at 834 G. Solid 
line: least square fit, arei oc a"^'^^"'*. Dashed line, theory: arei oc a^^-^^. The 
theoretical relaxation rate has been normalized to the experimental value at 
a = 78 nm. 
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cm^/s for the scattering length a « 110 nm. For K2 it is an order of magnitude 
higher at the same value of a [12] , which can be a consequence of a larger value 
of the characteristic radius of interaction R^. 

The suppression of the relaxation decay rate of weakly bound molecules of 
fermionic atoms has a crucial consequence for the physics of these molecules. 
At realistic temperatures the relaxation rate constant arei is much smaller than 
the rate constant of elastic collisions STra^^UT, where vt is the thermal velocity. 
For example, for the Li2 weakly bound molecules at a temperature T ~ 
and a ~ 8OOA, the corresponding ratio is of the order of 10"'* or 10"^. This 
opens wide possibilities for reaching BEC of the molecules and cooling the Bose- 
condensed gas to temperatures of the order of its chemical potential. Long-lived 
BEC of weakly bound molecules has been observed for ^°K2 at JILA [131 [HI 
and for ^Lia at Innsbruck [THUS], MIT [HI [18], ENS [19], Rice [20], Duke [21], 
and recently at Melbourne [63] and Tokyo [64]. Measurements of the molecule- 
molecule scattering length confirm the result add = 0.6a with the accuracy up 
to 30% [HdU. 

3 Heteronuclear molecules in Fermi-Fermi mix- 
tures 

3.1 Effect of mass ratio on elastic intermolecular interac- 
tion 

We now focus on novel physics of heteronuclear molecules which are expected to 
be formed in a mixture of two different fermionic atoms (Fermi-Fermi mixture) 
at a large positive 2-body scattering length a. In several aspects, the physics is 
similar to that discussed above for homonuclear molecules of fermionic atoms in 
different internal states. However, for a large mass ratio of the atoms, the situa- 
tion changes drastically. This is related to the existence of 3-body bound Efimov 
states, which in general makes it impossible to describe molecule-molecule scat- 
tering using only the value of the 2-body scattering length a. 

We start with calculating the amplitude of elastic interaction (scattering) 
between weakly bound heteronuclear molecules consisting of a heavy (mass M) 
and light (mass m) fermionic atoms, assuming that the atom-atom scattering 
length satisfies the inequality Re and again considering the ultracold limit 
determined by the condition (jl3p . In this case the scattering is dominated by 
the contribution of the s-wave channel, and we present here the exact results 
obtained in Ref. |65j using the zero-range approximation. Under the condition 
fca <C 1 the collision energy is much smaller than the molecular binding energy 
Eq. Hence, the s-wave molecule-molecule elastic scattering can be analyzed 
after setting the total energy equal to — 2£o = — fi.^//ia^. In the zero-range 
approximation one should solve the four-body free-particle Schrodinger equation 
which again can be written in the form (jlll) : 

- V^^ - V^, + * ^ 0, 
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Figure 7: Set of coordinates for the four-body problem with two heteronuclear 
molecules. 

where ri is the distance between two given heavy and light fermions, and r2 the 
distance between the other two (see Fig. 7). However, it is now more convenient 
to define the distance between the centers of mass of these pairs as /3R, and 
the separations between the heavy and light fermions in the other two possible 
heavy-light pairs as r± = a±ri + a;ipr2 ± /3R, with /3 = ^/^c^^aZ, = /^/-^^, 
a_ ~ fi/m, and fi = mM/(m + M) being the reduced mass. 
The symmetry condition (fTH]) then takes the form: 



*(ri, r2, R) = *(r2, ri, -R) = -*(r±, r^, ±/3(ri - rz) T («+ - a_)R), (34) 



and the Bethe-Peierls boundary condition should be applied for a vanishing 
distance in any pair of heavy and light fermions, i.e. for ri 0, r2 — » 0, and 
r-t ^ 0. For ri ^ it is again given by Eq. 

Due to the change in the definition of the coordinates the asymptotic ex- 
pression for the wavefunction ^E" at large distances R now reads: 



where the notation add is again used for the molecule-molecule scattering length, 
and the wavefunction of a weakly bound molecule is given by Eq. ([TS]). Then 
the asymptotic expression for the function /(r2,R) at large R is given by: 



For s-wave scattering the function / depends only on three variables: the 
absolute values of r2 and R, and the angle between them. Using the procedure 
described in Section 2 we obtain for / the same integral equation ([24| . The 
effect of different masses is contained in the expressions for the vectors S±^ 
which now read S± = {a=pr' ± /3R', a±r' =F /3R', =F/3r' =F (a+ - a_)R'}. 

In order to find the molecule-molecule scattering length as a function of the 
mass ratio M/m it is again more convenient to transform the integral equation 
for the function /(r, R) into an equation in the momentum space. Introducing 



^'«0o(ri)0o(r-2)(l-add//3i?); i? > a, 



(35) 



/(r2, R) w (2/r2a) exp (-r2/a)(l - add/f3R); i? > a. 



(36) 
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Figure 8: The ratio Udd/o. versus M/m. 



the Fourier transform of the function /(r, R) as /(k, p) = J (Prd^Rf{r, R) exp(zk • 
i(3p ■ R/a), we obtain the foUowing momentum-space equation: 



By making the substitution /(k, p) ~ {5{p) + g(k,p)/p^)/{l + k'^) we reduce 
Eq. (|37p to an inhomogeneous equation for the function g(k, p). This equation 
is similar to Eq. (pS)) and we do not present it here because of its complexity. 
As well as in the case of M = m, for p ^ the function ^(k, p) tends to a finite 
value independent of k. The molecule-molecule scattering length is again given 
by add = — 271^0 limp^o 3(k, p). In Fig. 7 we display the ratio add/a versus the 
mass ratio M/m, found in Ref. [65]. For the case of homonuclear molecules 
(m = M) we recover the molecule-molecule scattering length add = 0.6a. 

The universal dependence of add/a on the mass ratio, presented in Fig. 8, 
can be established in the zero-range approximation only if M/m is smaller than 
13.6. Calculations then show the absence of four-body weakly bound states, and 
for M/m ~ 1 the behavior of / suggests a soft-core repulsion between molecules, 
with a range ^ a. For the mass ratio larger than the limiting value 13.6, the 
description of the molecule-molecule scattering requires a three-body parameter 
coming from the short-range behavior of the three-body subsystem consisting 
of one light and two heavy fermions [53l [54j . A qualitative explanation of this 
behavior will be given in subsection 13.31 




/(k±a^(p^-p),pOdV 



2 +(k± a^(p'- p))2 -h (k± a±(p'-f- p))2 




/(k', -p) d3fc' 2^2(1 ^^2^ /3V)/(k, p) 



(37) 
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3.2 Collisional relaxation for moderate mass ratios 

The most exciting physics with weakly bound (heteronuclear) bosonic molecules 
consisting of different fermionic atoms is related to their collisional stability. As 
well as homonuclear molecules discussed above, they are in the highest rovibra- 
tional state and hence undergo relaxation into deeply bound states in molecule- 
molecule collisions, which leads to decay of the sample. The collisional relaxation 
determines the lifetime of the Bose gas of weakly bound molecules and there is a 
subtle question of whether and how the mass ratio M/ m can influence the sup- 
pression of this process [65] , originating from the Fermi statistics for the atoms 
and playing a crucial role in the case of homonuclear molecules. In a similar 
way, behaving themselves as point-like bosons at large intermolecular distances, 
heteronuclear molecules "remember" that they consist of fermions when the 
intermolecular separation becomes smaller than the molecule size (~ a). The 
relaxation requires the presence of at least three fermions at separations ~ i?e 
from each other. Two of them are necessarily identical, so that due to the Pauli 
exclusion principle the relaxation probability acquires a small factor propor- 
tional to a power of {qRe), where g 1/a is a characteristic momentum of the 
atoms in the weakly bound molecular state. What changes in this picture when 
the fermionic atoms forming the molecule have different masses? 

We first consider molecule-molecule relaxation collisions for the case where 
the mass ratio is smaller than the limiting value 13.6 and short-range physics is 
not supposed to influence the dependence of the relaxation rate on the 2-body 
scattering length a. As well as in the case of homonuclear molecules in Section 2, 
we assume the inequality a S> i?e and consider the ultracold limit described by 
Eq. (|13p. The conflguration space contributing to the relaxation probability 
can be again viewed as a system of only three atoms at short distances ^ 
Re from each other and a fourth atom separated from this system by a large 
distance ^ a. Hence, the four-body wavefunction decomposes into a product 
according to Eq. 5* = 77(2)^'^' (p, 17), with ^'(^^ being the wavefunction of 

the three-fermion system, and p <C a and fl being the hyperradius and the set of 
hyperangles for these fermions. The distance between their center of mass and 
the fourth atom is z, and the function 77(2) describes the motion of this atom. In 
the case of fermionic atoms with different masses there are two possible choices 
of a three-body subsystem out of four fermions. The most important is the 
relaxation in the system of one atom with the mass m and two heavier atoms 
with masses M. 

We then use the same arguments as in Section 2 and obtain equation ([5^ for 
the function at distances where i?e < p < a: ^'^^'> = A{a)^,y{fl)p''^^ , with 
the coefficient A{a) determining the o-dependence of the relaxation rate accord- 
ing to Eq. (I^T]) : arei — a'^' oc |y4(a)p. A similar scaling procedure as in Sec- 
tion 2 again leads to arei oc a"", where s — 21^+1, and restoring the dimensions 
we can write the relaxation rate in the form (j33p : Urei = C{hRe/rn){RelaY , 
with a coefficient C depending on the mass ratio and on short-range physics. 

However, the exponent s now depends not only on the symmetry of the 
three-body wave function but also on the mass ratio M/m. The smallest 
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Figure 9: The dependence of the exponent s = 2;^ + 1 in Eq. (j33|) on the mass 
ratio M/m. 

value of ly, i.e. the one corresponding to the leading relaxation channel at large 
a, is achieved for the p-wave symmetry in the system of one light and two heavy 
fermions [54|- In the interval — 1 < < 2 it is given by the root of the function 



where 7 = arcsin[M/(Af + m)]. A detailed derivation of Eq. ([55)1 is given in 
Ref. [Mj. 



In the case of equal masses we recover s = 2;/ + 1 w 2.55 obtained in 
Section 2, and it slowly decreases with increasing the mass ratio (see Fig. 9). 
For M/m ^ 1 nothing dramatic happens: the suppression of the relaxation 
rate with increasing the 2-body scattering length a becomes slightly weaker 
than for homonuclear molecules. However, for the mass ratio approaching the 
limiting value 13.6 the exponent s first reaches zero and then becomes negative, 
showing even an increase of the relaxation rate with a. We will give a qualitative 
explanation of this phenomenon using the Born-Oppenheimer approximation for 
the system of two heavy and one light atom. 

3.3 Born-Oppenheimer picture of collisional relaxation 

In the Born-Oppenheimer approximation one assumes that the state of a fast 
light atom adiabatically adjusts itself to the positions Ri and R2 of the slow 
heavy atoms. One then finds the wavcfunction and energy of a bound state 
of the light atom with two heavy atoms at a given separation between them, 
R — Ri — R2. For convenience, from this point on we change the notations and 




A(^) 



TTv i^sin7Cos(z^7 + 7) — sin(i/7) 



2 (z^ + 1) sin^ 7COS7sin(7ri^/2) 



(38) 
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use R to specify the distance between the heavy atoms, and r the coordinate 
of a Ught atom relative to their center of mass. In general, there are two states 
of a light atom in the field of two heavy ones: the gerade state (+) with the 
wavefunction remaining unchanged under permutation of the heavy atoms (R 
— R), and the ungerade state (— ) with the wavefunction changing its sign under 
this operation. The corresponding wavefunctions are given by 

- ( |r-R/2| ± |r + R/2| j ' ^''^ 

where J\f± are normalization coefficients which depend on R. The corresponding 
binding energies are 

e±(^) = -ft^K±(i?)/2m, (40) 
and the parameters k±(R) follow from the equation 

K±{R)Texp[-K±{R)R]/R=l/a. (41) 



Equation (j4ip is obtained by using the Bethe-Peierls boundary condition ^ 
for the wavefunctions of Eq. (|39p at vanishing light- heavy atom separations 
|r±R/2|. 

The ungerade (— ) state energy is always higher than the energy of the gerade 
(+) state. Moreover, for i? < a the ungerade state is no longer bound and we 
are dealing only with the gerade bound state. In the limit of i? <C a Eq. (j4ip 
gives K+ = 0.56. Then the energy of the gerade bound state, representing an 
effective potential for the relative motion of the heavy atoms, is given by 

e+(i?) = -OAeh^/mR^. (42) 

We thus see that when the heavy atoms are separated from each other by a 
distance R ^ a, the light atom mediates an effective 1/R^ attraction between 
them. Actually, the same result follows from the Efimov picture of effective 
interaction in a three-body system 53J and the Born-Oppenheimer approxima- 
tion only gives a physically transparent illustration of this picture |651 166] . For 
a large mass ratio the mediated attractive potential e+(i?) = — 0.16?i^/mi?^ 
strongly modifies the physics of the relaxation process. It competes with the 
Pauli principle which in terms of effective interaction manifests itself in the cen- 
trifugal 1/ Rl^ repulsion between the heavy atoms. The presence of this repulsion 
is clearly seen from the fact that the light-atom wavefunction V'r(i") does not 
change sign under permutation of heavy fermions. As the total wavefunction 
of the three-body system -(/'^(r)x(R) is antisymmetric with respect to this per- 
mutation, the wavefunction of the relative motion of heavy atoms x(R) should 
change its sign. Therefore, x(I^) contains only partial waves with odd angular 
momenta, and for the lowest angular momentum (p-wave) the centrifugal barrier 
is Uc(R) = 2h^/MR^. For comparable masses it is significantly stronger than 
e+(i?). Thus, we have the physical picture discussed in the case of homonuclear 
molecules: the Pauli principle (centrifugal barrier) reduces the probability for 
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the atoms to be at short distances and, as a consequence, the relaxation rate 
decreases with increasing the atom-atom scattering length a. 

The role of the effective attraction increases with M/ra. As a result, the 
decrease of the relaxation rate with increasing a becomes weaker. The exponent 
s in Eq. (|33p continuously decreases with increasing M/m and becomes zero 
for M/m = 12.33 (see Fig. 9). In the Born-Oppenheimer picture this means 
that at this point one has a balance between the mediated attraction and the 
centrifugal repulsion. A further increase in M/m makes s negative and it reaches 
the value s = — 1 for the critical mass ratio M/m — 13.6. Thus, in the range 
12.33 < M/m < 13.6 the relaxation rate increases with a. 

For an overcritical mass ratio M/m > 13.6 we have a well-known phe- 
nomenon of the fall of a particle to the center in an attractive 1/R^ potential 
[25] . In this case the shape of the wavefunction at distances of the order of Re 
can significantly influence the large-scale behavior and a short-range three-body 
parameter is required to describe the system. The wavefunction of heavy atoms 
x(R) acquires many nodes at short distances R, which indicates the appearance 
of 3-body bound Efimov states. 

3.4 Molecules of heavy and light fermionic atoms 

The discussion of the previous subsection shows that weakly bound molecules of 
heavy and light ferniions become coUisionally unstable for the mass ratio M /m 
close to the limiting value 13.6. The effect of the Pauli principle becomes weaker 
than the attraction between heavy atoms at distances R a, mediated by light 
fermions. However, this picture explains only the dependence of the relaxation 
rate on the 2-body scattering length a. At the same time, for heteronuclear 
molecules the relaxation rate and the amplitude of elastic molecule-molecule 
interaction can also depend on the mass ratio irrespective of the value of a and 
short-range physics. To elucidate this dependence, we will look at the interac- 
tion between the molecules of heavy and light fermions at large intermolecular 
separations. 

We consider the interaction between two such molecules in the Born-Oppenheimer 
approximation and calculate the wavefunctions and binding energies of two light 
fermions in the field of two heavy atoms fixed at their positions Ri and R2. 
The sum of the corresponding binding energies gives an effective interaction po- 
tential Ucs for the heavy fermions as a function of the separation R = |Ri — R2I 
between them. 

For R > a, there are two bound states, gerade {+) and ungerade (-), for a 
light atom interacting with a pair of fixed heavy atoms. Their wavefunctions are 
given by Eq. ([55]) . and the corresponding binding energies follow from Eqs. (|40p 
and ([iT|). For large R satisfying the condition exp(— i?/a) <C 1, Eq. (|^ yields: 

e±{R} « ~\eo\ T 2|eo|| exp(-i?/a) + (43) 
where the binding energy of a single molecule, eg, is given by Eq. ^ with the 
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Figure 10: Interaction potential for two molecules of heavy and light fermions 
as a function of the separation R between the heavy atoms. 



reduced mass very close to the light atom mass m. 

UUR) - 4|6o|| (l - ^) exp(-2J?/a). (44) 

Since the light fermions are identical, their two-body wavefunction is an 
antisymmetrized product of gerade and ungerade wavefunctions; 

V'R(ri,r2) = [V^+CrOV-RCra) - V^^(r2)VR(ri)]/\/2. (45) 

The Born-Oppcnheimer adiabatic approach is valid at distances R > a, where 
the effective interaction potential between the molecules, Ucs, is the sum of 
e+(i?) and e„(i?) (one should add 2|eo| so that Uos{R) — > for R — > oo). This 
potential is displayed in Fig 10. 

At sufficiently large inter-heavy separations where Eq. is valid, the ef- 
fective potential can be written as 

f/eff(i?) = e+{R) + e_(i?) + 2|eo| « U,^{R). (46) 

The potential Ucx originates from the exchange of light fermions and thus can 
be treated as an exchange interaction. It is purely repulsive and, according to 
Eq. (UH), has the asymptotic shape of a Yukawa potential at large R. Direct 
calculations show that Ucx is a very good approximation to [/off for R > 1.5a. 

We now demonstrate the calculation of the dimer-dimer scattering length 
add in the limit of M/m ^ 1 |67] . In the Born-Oppenheimer approach the 
Schrodinger equation for the relative motion of two molecules reads; 

(-(?iV"^)V^ + C/off(i?)-e)*(R) -0, (47) 

where e is the collision energy. Note, that the repulsive effective potential is in- 
versely proportional to the light mass m, whereas the kinetic energy operator in 
Eq. (IT71) has a prefactor 1/M. Therefore, for a large mass ratio M/m, the heavy 
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atoms approach each other at distances sraaller than a with an exponentially 
small tunneling probability P cx exp(— B-y/M/m), where B ^ 1. This leads to 
the relaxation rate constant 



arei oc exp (- B yjM/m) , (48) 

which strongly decreases with increasing the mass ratio M/m. 

The analysis shows that the elastic part of the scattering amplitude can 
be calculated with a very high accuracy using Eq. (|T7)) for M/m > 20 and 
is practically insensitive to the way we choose the boundary condition for the 
wavefunction at i? = a. The dominant contribution to the scattering comes 
from distances in the vicinity of R — add 3> a, where the effective potential can 
be approximated by Eq. ()44|) with a constant preexponential factor: 



f/eff (i?) « 2h^{maaMr^ exp(-2i?/a). (49) 
Then, the zero energy solution of Eq. (|47l) that decays at smaller R, reads: 



*(i?) = >K/--e-/M, (50) 
R y\ m add J 

where Kq is the decaying Bessel function. Comparing the result of Eq. ([50)) at 
large R with the asymptotic behavior ^'(i?) (x (l — add/R) we obtain an equation 
for add- 

a_^Je^Mj^\ 



(^dd^^\n\- . (51) 

2 \ 2 m addj 

This gives 

add w a In i/M/m, (52) 

and the scattering cross section is 

add = Stto^^. (53) 

From Eq. ([50]) we see that the interval of distances near R — add, where the 
wavefunction changes, is of the order of a. This justifies the use of Eq. 
In fact, the corrections to Eq. ([STjl can be obtained by treating the difference 
between Eqs. ([44]) and ([49]) perturbatively. In this way the first order correction 
to the dimer-dimer scattering length is —{3/4:)a'^/add, where add is determined 
from Eq. ([5T|) . 

Qualitatively, Ucs (R) can be viewed as a hard-core potential with the radius 
add, where the edge is smeared out on a lengthscale ~ a <C add- Therefore, the 
ultracold limit for dimer-dimer collisions, required for the validity of Eq. ([53)) . 
is realized for relative momenta of the dimers, k, satisfying the inequality 

kadd < 1. (54) 

It can be useful (see Ref. [68]) to approximate the potential UeS by a pure hard 
core with the radius add- This approximation works under the condition /ca <C 1, 
which is less strict than Eq. ([M]) . 
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Let us now mention that numerical calculations Ref. [S7] show that there are 
no resonances in the dimer-dimer scattering amplitude, which could appear in 
the presence of a weakly bound state of two dimers. Here we give a qualitative 
explanation of the absence of these bound states. Suppose there is such a state 
with energy e — > 0. Then, at distances R > a the wavefunction of the heavy 
atoms should exponentially decay on the distance scale ~ a^mjM <C a, since 
C/eff represents a barrier with the height ^ 1/ma?. This means that the heavy 
atoms in such a bound state should be localized mostly at distances smaller 
than a. The gerade light atom is also localized at these distances as seen from 
the shape of the function -0"*". The motion of the ungerade light atom relative to 
the localized trimer can be viewed as scattering with odd values of the angular 
momentum, and due to the centrifugal barrier the bound states of this atom 
with the trimer should be localized at distances ^ a from the heavy atoms. In 
this case one would expect the BO approximation to work, since the ungerade 
light atom is moving much faster than the heavy atoms. However, this leads to a 
contradiction, because in the BO approach discussed above the ungerade state 
at interheavy separations i? < a is unbound. We thus conclude that weakly 
bound states of two dimers are absent. 

Although there are no resonances in the dimer-dimer collisions, there are 
branch-cut singularities in the scattering amplitude. They are related to the 
presence of inelastic processes in molecule-molecule collisions. These represent 
the relaxation of one of the colliding dimers into a deeply bound state, the 
other dimer being dissociated, and the formation of bound trimers consisting of 
two heavy and one light atom, the other light atom carrying away the released 
binding energy. 

3.5 Trimer states 

The trimer states, which in most cases can be called Efimov trimers, are in- 
teresting objects. Their existence can be seen from the BO picture for two 
heavy atoms and one light atom in the gerade state. Within the BO approach 
the three-body problem reduces to the calculation of the relative motion of the 
heavy atoms in the effective potential created by the light atom. For the light 
atom in the gerade state, this potential is e+(i?) found in the previous subsec- 
tion. The Schrodinger equation for the wavefunction of the relative motion of 
the heavy atoms, Xy(R), reads: 

i?x(R) = [-{re/M)Vl, + e+(i?)] x.(R) = e,x.(R). (55) 

The trimer states are nothing else than the bound states of heavy atoms in the 
effective potential e+(i?). Accordingly, they correspond to the discrete part of 
the spectrum Cj^, where the symbol v denotes a set containing angular {I) and 
radial (n) quantum numbers. For R <^ a the potential e+(i?) is proportional to 
— (see Eq. (|42p ) and, if this effective attraction overcomes the centrifugal 
barrier, we arrive at the well known phenomenon of the fall of a particle to 
the center in an attractive 1/R^ potential. Then, for a given orbital angular 
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momentum Z, the radial part of Xf can be written as 



Xy{R)<^R ^/^sin(s/lni?/ro), R 'ti a 



(56) 



where 



si = v^O.lGAf/m- {I + 1/2)2. 



(57) 



The three-body parameter vq determines the phase of the wavefunction at small 
distances and, in principle, depends on /. The wavefunction (|56p has infinitely 
many nodes, which means that in the zero-range approximation there are in- 
finitely many trimer states. This is one of the properties of three-body systems 
with resonant interactions discovered by Efimov [53]. We see that the fall to 
the center is possible in many angular momentum channels, provided the mass 
ratio is sufficiently large. However, for practical purposes and for simplicity, 
it is sufficient to consider the case where the Efimov effect occurs only for the 
angular momentum channel with the lowest possible I for a given symmetry. 
This implies that when the heavy atoms are fermions and one has odd I, in 
order to confine ourselves to Z = 1 we should have the mass ratio in the range 
14 < M/ni < 76. For bosonic heavy atoms where / is even, we set I = and 
consider M/m < 39 to avoid the Efimov effect for I > 2. In both cases we need 
a single three-body parameter tq. 

The formation of Efimov trimers in ultracold dimer-dimer collisions is ener- 
getically allowed only if < — 2|eo|. This means that the trimers that we are 
interested in are relatively well bound and their size is smaller than a. There- 
fore, the process of the trimer formation is exponentially reduced for large mass 
ratios as the heavy atoms have to tunnel under the repulsive barrier [/off(R). 
Moreover, this process requires all of the four atoms to approach each other at 
distances smaller than a, and its rate decreases with the trimer size because it 
is more difficult for two identical light fermions to be in a smaller volume. 

From Eq. one sees that the behavior of the three-body system does not 
change if rg is multiplied by 



On the other hand, the dimensional analysis shows that the quantity ejy/eo 
depends only on the ratio a/ro. This means that except for a straightforward 
scaling with a, properties of the three-body system do not change when a is mul- 
tiplied or divided by A/ . This discrete scaling symmetry of a three-body system, 
which shows itself in the log-periodic dependence of three-body observables, has 
yet to be observed experimentally. In the case of three identical bosons, where 
the BO approach does not work and it is necessary to solve the three-body 
problem exactly [53^ , the observation of the consequences of the discrete scaling 
requires to change a by a factor of A w 22.7, which is technically very difficult in 
ongoing experiments with cold atoms. In this respect three-body systems with 
a very large mass difference can be more favorable because of smaller values of 
A. For example, in order to see one period of the log-periodic dependence in a 
Cs-Cs-Li three-body system a has to be changed only by a factor of A « 5. 



A/ = exp(7r/s;). 



(58) 
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At this point it is worth emphasizing that three-body effects can be observed 
in a gas of Ught-heavy dimers, where the interdimer repulsion originating from 
the exchange of the hght fermions strongly reduces the decay rate associated 
with relaxation of the dimers into deep bound states. The trimer formation 
in dimer-dimer collisions is very sensitive to the positions and sizes of Efimov 
states, and the measurement of the formation rate can be used to demonstrate 
the discrete scaling symmetry of a three-body system. Indeed, this rate should 
have the log-periodic dependence on a and is detectable by measuring the life- 
time of the gas of dimers. 

Besides the Efimov trimers, one can have of one light and two heavy atoms 
may form "universal" trimer states well described in the zero-range approxi- 
mation without introducing the three-body parameter [69]. In particular, they 
exist for the orbital angular momentum I — 1 and mass ratios below the critical 
value, where the Efimov effect is absent and short-range physics drops out of 
consideration. One of such states emerges at M/m « 8 and crosses the trimer 
formation threshold {etr = ~2|eo|) at M/m « 12.7. The existence of this state 
is already seen in the BO picture. It appears as a bound state of fermionic heavy 
atoms in the potential e+(i?) for ^ = 1. The other state exists at M/m even 
closer to the critical mass ratio and never becomes sufficiently deeply bound to 
be formed in cold dimer-dimer collisions. The universal trimer states also exist 
for I > 1 and M/m > 13.6 However, the trimer formation in dimer-dimer 
collisions at such mass ratios is dominated by the contribution of Efimov trimers 
with smaller I. Therefore, below we focus on the formation of Efimov trimers. 

The calculation of the intrinsic lifetime of a trimer requires a detailed knowl- 
edge of short-range physics and is a tedious task. Estimates of the imaginary 
part of the trimer energy, r~^, from the experimental data on Csa trimers [70] 
show that it is approximately by a factor of 4 smaller than the real part (in 
this case ry* « 0.06). From a general point of view, we do not expect that the 
trimers with a binding energy Ci, < — 2|eo| are very long-lived. However, one can 
have relatively narrow resonances, and we demonstrate calculations for various 
values of the elasticity parameter 77* [67] . 

3.6 Collisional relaxation of molecules of heavy and light 
fermions and formation of trimers 

Let us now discuss inelastic processes in dimer-dimer collisions and start with 
relaxation of the dimers into deeply bound states. The typical size of a deeply 
bound state is of the order of the characteristic radius of the corresponding 
interatomic potential. We first consider the relaxation channel that requires one 
light and two heavy atoms to approach each other at distances a. Unlike 

the trimer formation, this decay mechanism is a purely three-body process. The 
other light atom is just a spectator. A qualitative scenario of this process is the 
following. With the tunneling probability which is exponentially suppressed for 
large M/m, two dimers approach each other at distances R ^ a. Then the 
heavy atoms are accelerated towards each other in the potential e+(i?), and the 
light atom in the gerade state is always closely bound to the heavy ones as is 
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seen from the shape of the function . The relaxation transition occurs when 
the heavy atoms (and the gerade Hght fermion) are at interatomic separations 
~ i?e- The relaxation rate constant certainly satisfies Eq. (|48p . but we should 
also find out how to take into account the relaxation process in the description 
of the trimer states and their formation. The most convenient way to do so is 
to consider the three-body parameter rg as a complex quantity and introduce 
the so-called elasticity parameter 77* — — s;Arg(ro) [71]. As follows from the 
asymptotic expression for the wavefunction a negative argument of ro 

ensures that the incoming flux of heavy atoms is not smaller than the outgoing 
one: 

*out/$i„ = exp[4s,Arg(ro)] = cxp[-4?7*] < 1. (59) 

This mimics the loss of atoms at small distances due to the relaxation into 
deeply bound states. In the analysis of Efimov states, the imaginary part of 
ro leads to the appearance of an imaginary part of tv. This means that any 
Efimov state has a finite lifetime r due to the relaxation. For small |Arg(ro)| 
and for trimer states that are localized at distances smaller than a, we get 
/\ev\ = 4|Arg(ro)| = Arj^/si. Strictly speaking this fact indicates that it 
is not possible to separate the relaxation process from the trimer formation 
because the trimers that are formed in dimer-dimer collisions will eventually 
decay due to the relaxation. Nevertheless, both the modulus and the argument 
of the three-body parameter can be determined by measuring the lifetime of a 
gas of dimers, leading to a number of quantitative predictions concerning the 
structure of Efimov states in the three-body subsystem of one light and two 
heavy atoms. 

Another relaxation channel is the one in which two light atoms approach a 
heavy atom at distances Re <^ a. This channel is, however, suppressed due to 
the Fermi statistics for the light atoms, which strongly reduces the probability 
of having them in a small volume. As a result, for realistic parameters this 
relaxation mechanism is much weaker than the one in the system of one light 
and two heavy atoms [67J. 

The study of the formation of trimer states in molecule-molecule collisions 
requires us to go beyond the conventional Born-Oppenheimer approximation, 
since this approximation breaks down for the ungerade light atom at separa- 
tions between the heavy atoms R < a. Within the recently developed "hybrid 
Born-Oppenheimer" approach j67j one applies the Born-Oppenheimer method 
to the gerade light fermion which is characterized by the wavefunction ■0R,(r) 
and energy e+(i?) adiabatically ajusting themselves to the motion of the heavy 
atoms. The gerade light atom is then integrated out by introducing the poten- 
tial €+{R) for the heavy atoms. Once this is done the original 4-body problem 
is reduced to a 3-body problem described by the Schrodinger equation equation 

[H~{ny2fi3)Wl-E]^{R,r) = 0, (60) 

where H is given by Eq. dSS]), ^3 = 2mM/{2M + m), E = -2|eo| +e is the total 
energy of the four-body system in the center of mass reference frame, and e is 
the dimer-dimer collision energy. This problem is then treated exactly in the 
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Bethe-Peierls approach, and the interaction of the hght atom with the heavy 
atoms is included in the form of the Bethe-Peierls boundary condition ([8]) for 
^ at vanishing light-heavy separations |r ± R/2|. The ungerade symmetry for 
this atom is taken into account by the condition 

*(R,r) = -^(R.-r). (61) 

As the heavy atoms are identical fermions, we have 4'(R, r) = — R, r). 
Combined with Eq. ([6T|) . this leads to the condition ^'(R, r) — ^'(— R, — r). 
Therefore, ^(R, r) describes atom-dimer scattering with even angular momenta, 
and for ultracold collisions we have to solve an s-wave atom-dimer scattering 
problem. 

In order to solve Eq. (|60|) we follow the method of Ref. [54 and introduce an 
auxiliary function /(R) and write down the wavefunction ^'(R, r) in the form: 

vI/(R,r)-V/ x.(R)x:(R')^«.(2r,R')/(R'), (62) 



where 



and 



g-K„|r-R72| g-K„|r+R72| 

^-(^-'^^)^ 4.|r-R/2r4.|r + R72| ^''^ 



(64) 



" \ ^i^2ii;iE^7;;)/h, < e. 

For < E, the trimer can be formed in the state ly. In this case k^, is 
imaginary and the function (|63p describes an outgoing wave of the light atom 
moving away from the trimer. The choice of the sign in Eq. ([M)) ensures that 
there is no incoming flux in the atom-trimer channel. 

Using the Bethe-Peierls boundary condition ([8]) for the wavefunction (|62p at 
|r ± R/2| we one obtain an integral equation for the function /(R): 



^2J2^i{eo-E)/h-l/a 



L-L' + sin^ g ^ " ^ ^ /(R) = 0, (65) 



where fi = mM/{m + M), 9 = arctan ^1 + 2M/m, and 

L'f iR)=[ ^[x.(R)x:(R')i^K.(R,R')-x"(R)x°*(R')^.o(R,R')]/(R'), (66) 

LfiK)=pj |g(|R-R'|) [/(R)-/(R')]±G(v/i?2 + i?'2_2R.Rcos26i)/(R')}, (67) 

sin2gAf(£o - E)K2WM[e^ - E)X/hsm9) 
GiX) , (68) 

with K2{z) being the exponentially decaying Bessel (Macdonald) function. A 
detailed derivation is given in Ref. [S7] and is omitted here. The operators L and 
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L' conserve angular momentum and, expanding the function /(R) in spherical 
harmonics, we arrive at a set of uncoupled ID integral equations for each of 
the radial functions fi{R). Below we present the results for s-wave dimer-dimer 
scattering. 

At large distances {R ^ a) the reduced wavefunction ^'(R, r) takes the form: 
vI/(R,r)«vl/(R)^^(r), (69) 

and one can show that 

^'(R)cx/(R); i?>a. (70) 

Therefore, /(R) can serve as the wavefunction for the dimer-dimer motion at 
large distances. In particular, it contains the dimer-dimer scattering phase shift. 

The dimer-dimer s-wave scattering amplitude add is determined from the 
asymptotic behavior of the solution of Eq. ([65)1 at large distances for E = 2eo, 
which should be matched with 

fo{R)<^{l/R-l/add) (71) 

at i? 3> a\n ^yM/rn. In Fig. 11 we compare the resulting add/o, with that 
following from Eq. ([5T|) . The results agree quite well even for moderate values 
of M/m. These results also agree with the calculations based on the exact four- 
body equation for M/m < 13.6 |65j and displayed in Fig. 8, and with the Monte 
Carlo results for M/m < 20 [72]. 
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Figure 11: The dimer-dimer s-wave scattering length add/ a. The solid curve 
shows the results obtained in the hybrid Born-Oppenheimer approximation, and 
the dotted line the results of Eq. ((5T|) . 

It is straightforward to extend this theory to account for inelastic processes 
of the trimer formation and the relaxation of dimers into deeply bound states. 
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Let us first assume that the rate of the relaxation into deep molecular states is 
negligible and neglect this process. Then the three-body parameter is real, and 
the trimer formation rate is determined by the imaginary part of the s-wave 
scattering length. The rate constant is given by [5S] 

16hTT 

a = ^lm(add). (72) 

Alternatively, if it is necessary to know the rate of the trimer formation in the 
state ly, one can substitute the solution of Eq. (|65|) into Eq. (|62p and calculate 
the flux of light atoms at r — s- cx). The summation over v gives the same result 
as Eq. ([7^ . We find that the contribution of the highest "dangerous" trimer 
state is by far dominant and a is very sensitive to its position. 

We now include the relaxation of the dimers into deeply bound states. As we 
have mentioned above, the light-heavy-heavy relaxation process can be taken 
into account by adding an imaginary part to the three-body parameter. The to- 
tal inelastic decay rate is then still given by Eq. (|7^ . However, strictly speaking 
we can no longer distinguish between the formation of the trimer in a particular 
state and the collisional relaxation since the trimers ultimately decay due to 
the relaxation process. In this sense the only decay channel is the relaxation. 
However, for a sufficiently long lifetime of a trimer, i.e. if the trimer states are 
narrow resonances, we can still observe a pronounced dependence of the total 
inelastic decay rate on the position of the highest "dangerous" trimer state (see 
below). 

Fig. [T^ shows the results for the inelastic collisional rate in the case of bosonic 
molecules with the mass ratio M/m — 28.5 characteristic of ^^^Yb-^Li dimers. 
The solid line corresponds to the case of a real three-body parameter. It is 
convenient to introduce a related quantity, uq, defined as the value of a at 
which the energy, e^, of a trimer state exactly equals E — 2eo. This new 
"dangerous" trimer state becomes more deeply bound for a > oq and the rate 
constant rapidly increases. It is proportional to the density of states in the 
outgoing atom-trimer channel. The corresponding orbital angular momentum 
is equal to 1 and the threshold law reads (see also inset in Fig. [T2|l : 

a (X const + {E - tyf^"^ oc const + (a - aafl"^. (73) 

The constant term in Eq. (j73|) describes the contribution of more deeply bound 
states, which is typically very small. In fact, as trimer states become more 
compact, both light atoms should approach the heavy atoms and each other to 
small distances where the trimer formation takes place. Since they are identical 
fermions, there is a strong suppression of the trimer formation to these deeply 
bound states. 

The dependence of a on a/ao is periodic on the logarithmic scale, the mul- 
tiplicative factor being equal to Ai « 7.3. The dashed, dotted, and dash-dotted 
curves are obtained for 77, =0.1, 0.5, and 1, respectively. The corresponding val- 
ues of the ratio $out/*i'in are 0.67, 0.14, and 0.02. The horizontal line represents 
the limiting case of 77, = 00 or $out = 0. This case is universal in the sense 
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Figure 12: The inelastic rate constant for bosonic dimers with M/m = 28.5 as 
a function of the atom-atom scattering length a. The solid line corresponds to 
the case of a real three-body parameter. The results plotted in dashed, dotted, 
dash-dotted, and dash-dot-dot lines are obtained by taking into account the 
light-heavy-heavy relaxation processes. The values of the elasticity parameter 
77* =0.1, 0.5, 1, and cxd, respectively (see text). The quantity oq is the value 
of a at which the energy of a trimer state equals E = 2eo and a new inelastic 
channel opens. The inset shows the region o w ao in greater detail in order to 
see the threshold behavior ([75)1 . 



that physical observables depend only on the masses and the atomic scattering 
length. 

For a very weak light-heavy-heavy relaxation, the dimer-dimer inelastic col- 
lision can be viewed as the formation of a trimer (with the rate constant a) 
followed by its slow decay due to the relaxation. In this case one can think 
of detecting the trimers spectroscopically. We note, however, that even for the 
conditions corresponding to the dashed curve in Fig. 1121 i.e. for 77* as small as 
0.1, the decay rate of the trimer t~'^ k, 0.25|e^|//i is rather fast, which will likely 
make its direct detection difficult. For larger ry* it is impossible to separate the 
formation of trimers from their intrinsic relaxational decay, and a is practically 
the relaxation rate constant. Remarkably, it remains quite sensitive to the po- 
sitions of trimer states (in this case resonances) for values of 77* up to 0.5 and 
even larger. This suggests that measuring the lifetime of a gas of dimers as 
a function of a may provide important information on three-body observables. 
Moreover, for small 77* it may be possible to crate a stable molecular gas in 
sufficiently broad regions of a, where "dangerous" trimer states are far from the 
trimer formation threshold. 
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The inelastic rate a for other mass ratios in the range 20 < M/m < 76 has 
also been found [67] . Its dependence on the scattering length has the same form 
as depicted in Fig. [T^l The maximum of the rate constant is well fitted by the 
formula amax = 5.8{ha/M) exp(— 0.87-\/M/m), and the position of the horizon- 
tal line (77* = 00) by aoo = l-6{ha/M) exp(— 0.82-\/Af/m). The multiplicative 
factor in the log-periodic dependence is given by Eqs. ((58)) and (|57|) with ^ = 1. 

The same method was employed to estimate the formation of the universal 
trimer state with the orbital angular momentum Z = 1 at mass ratios M/m > 
12.7 but below the critical value for the onset of the Efimov effect [ST]. The rate 
constant increases with M/m and reaches a — 0.2{ha/M) close to the critical 
mass ratio. This corresponds to the imaginary part of the scattering length 
Imadd ~ 4 X 10~^a, which is by a factor of 300 smaller than the real part of 
add obtained from four-body calculations [65] . Thus, the formation of this state 
does not change the elastic scattering amplitude add shown in Fig. 8. 

We can now estimate the coUisional rates for ^^^Yb-^Li dimers. On the 
basis of the results in Fig. [H] we find that for a = 20nm the upper bound of the 
inelastic rate constant is oimax ~ 4 x 10~^^cm'^/s. The elastic rate constant for 
a thermal gas with a ~ 20nm and T ^ lOOnK equals 

aei « 8n\add\W^T/M ~ 4 x 10-"cmVs. (74) 

Here we used the calculated s-wave dimer-dimer scattering length add ~ 1.4a. 
We see that aei is much larger than a and this inequality becomes even more 
pronounced for larger a due to the scaling relations aei oa and amax oc 
a. Thus, the gas of molecules of heavy and light fermions is well suited for 
evaporative cooling towards their Bose-Einstein condensation. 

4 Crystalline molecular phase 

4.1 Born-Oppenheimer potential in a many-body system 
of molecules of heavy and light fermions 

Strong long-distance repulsive interaction between weakly bound molecules of 
light and heavy fermionic atoms has important consequence not only for the 
relaxation process, but also for macroscopic properties of the molecular system. 
In contrast to two-component Fermi gases of atoms in different internal states, 
heteronuclear Fermi- Fermi mixtures can form a molecular crystalline phase even 
when the mean interparticle separation greatly exceeds the size of the molecule: 

i? > a. (75) 

Let us consider a mixture of heavy and light fermionic atoms with equal 
concentrations and a large positive scattering length for the interaction be- 
tween them, satisfying the inequality a ^ Re- At zero temperature all atoms 
will be converted into weakly bound molecules and under the condition (|75p 
the molecular size {'^ a) will be much smaller than the mean intermolecular 
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separation. Then, using the Born-Oppenheimer approximation and integrating 
out the motion of light atoms we are left with a system of identical (composite) 
bosons which is described by the Hamiltonian: 

^ = E + ^ E UM,), (76) 

i i^j 

where indices i and j label the bosons, their coordinates are denoted by and 
Rj, and Rij = |Ri — Rj| is the separation between the i-th and j-th bosons. 
Assuming that the motion of light fermions is three-dimensional, the effective 
repulsive potential UcS is given by Eq. (|44)) and is independent of the mass of 
the heavy atom M. Therefore, at a large mass ratio M/m it dominates over 
the kinetic energy which is inversly proportional to M, which may lead to the 
formation of a crystalline phase. 

We will discuss the case where the motion of heavy atoms is confined to 
two dimensions, while the motion of light atoms can be either 2D or 3D. It 
will be shown that the Hamiltonian ([76]) with U^s (El) , supports the first order 
quantum gas-crystal transition at T = [68j . This phase transition resembles 
the one for the flux lattice melting in superconductors, where the flux lines are 
mapped onto a system of bosons interacting via a 2D Yukawa potential j73| . 
In this case the Monte Carlo studies [TH [7S] identified the first order liquid- 
crystal transition at zero and finite temperatures. Aside from the difference in 
the interaction potentials, a distinguished feature of our system is related to its 
stability. The molecules can undergo coUisional relaxation into deeply bound 
states, or form weakly bound trimers. Another subtle question is how dilute 
the system should be to enable the use of the binary approximation for the 
molecule-molecule interaction, leading to Eqs. (|76|) and (|44)l . 

Let us first consider the system of N molecules and derive the Born-Oppenheimer 
interaction potential for this system. Omitting the interaction between light 
(identical) fermions, it is sufficient to find N lowest single-particle eigenstates, 
and the sum of their energies will give the interaction potential for the molecules. 
For the interaction between light and heavy atoms we use the Bethe-Peierls ap- 
proach, and the wavefunction of a single light atom then reads: 

N 

*({R},r) = ^C,G,(r-R,), (77) 

i=i 

where r is its coordinate. The Green function satisfies the equation (— -I- 
k^)Gk,{y) — 5{y). The energy of the state ((77)) equals e = — /i^K^/2m, and 
here we only search for negative single-particle energies. The dependence of 
the coefficients d and k on {R} is obtained using the Bethe-Peierls boundary 
condition: 

^'({R},r)cxG,„(r-R,); r ^ R,. (78) 

Up to a normalization constant, G^o is the wavefunction of a bound state of 
a single molecule with energy eo = —h? K'^/2m and molecular size Kq^. The 
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Bethe-Peierls boundary condition ([75| written through the Green function G^o 
can be used for both 2D and 3D motion of Ught atoms. In the latter case one 
has kq = a~^, the Green function is GKo(r ~ R-i) oc (|r — Ri|^^ — a~^) for 
|r - R,| ^ 0, and Eq. ^ takes the form ®. 

From Eqs. ([77)1 and ((75|) we get a set of N equations: J^j ^ij^j = 0' where 
Aij = A(K)(5y + GK(i?ij)(l - (5,y), = |Ri - Rjl, and A(k) = \m\r^Q[G^(r) - 
Gko(^)]- The single-particle energy levels are determined by the equation 

det[Ay(K,{R})] = 0. (79) 

For Rij oo, Eq. ([79]) gives an N-fold degenerate ground state with k ^ kq. 
At finite large Rij, the levels split into a narrow band. Given a small parameter 

^ = G«„(i?)/Ato|A'jKo)| «1, (80) 

where i? is a characteristic distance at which heavy atoms can approach each 
other, the bandwidth is Ae w 4|eo|C ^ ko|- It is important for the adiabatic 
approximation that all lowest N eigenstates have negative energies and are 
separated from the continuum by a gap ~ |eo|. 

We now calculate the single-particle energies up to second order in ^. To 
this order we write k(A) « kq + k'^A -I- k'^^X^ /2 and turn from Aij(K) to Aij(X): 

A^j = X6,j + [G^„ (R,,) + K\XdG^„ {R,,)/dK\ (1 - 5,,), (81) 

where all derivatives are taken at A = 0. Using Aij ([5T|) in Eq. ([75)1 gives a 
polynomial of degree N in A. Its roots A^ give the light-atom energy spectrum 
ei = — ?i^K^(A,;)/2m. The total energy, E = Ylif=i then given by 

N N 

1=1 1=1 

Keeping only the terms up to second order in ^ and using basic properties of 
determinants and polynomial roots we find that the first order terms vanish, 
and the energy reads E ~ Neo + (1/2) J2i^j where 

U{R) = -^L U«)2^^ + («<)',G^^(i?) . (83) 

Thus, up to second order in ^ the interaction in the system of N molecules is 
the sum of binary potentials ([83]) . 

If the motion of light atoms is 3D, the Green function is Gk{R) = {I/AttR) exp(- 
and A(k) = (kq ~ '«)/47r, with the molecular size Kq ^ equal to the 3D scattering 
length a. Equation then gives a repulsive potential ?7ex (|^ which we now 
denote as U^^: 

Ul^[R) - 4|eo|(l - (2«;oi?)-')exp(-2Koi?)Aoi?, (84) 

and the criterion ((80|) reads (1/Koi?) exp(— Koi?) <C 1. For the 2D motion of 
light atoms we have Gk,{R) = (l/27r)isro(Ki?) and X(k) = — (l/27r) ln(K/Ko), 
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where Kq is the decaying Bessel function, and Kq ^ follows from [73] . This leads 
to a repulsive intermolecular potential: 

C/2f (i?) = 4\eo\[KoRKoiKoR)Ki{KoR) - Ki{KaR)], (85) 

with the validity criterion Ko{koR) ^ 1. In both cases, which we denote 2x3 
and 2x2 for brevity, the validity criteria are well satisfied already for kqR « 2. 

4.2 Gas-crystal quantum transition 

The inequality kqR > 2 may be considered as the condition under which the 
system is described by the Hamiltonian (|76p. with C/e// given by Eq. ((84| for 
the 3D motion of light atoms, or by Eq. ([55)1 in the case where this motion 
is 2D. The state of the system is then determined by two parameters: the 
mass ratio M/m and the rescaled 2D density ukq^. At a large M/m, the 
potential repulsion dominates over the kinetic energy, which should lead to 
the formation of a crystalline ground state. For separations Rij < Kq^ the 
adiabatic approximation breaks down. However, the interaction potential U{R) 
is strongly repulsive at larger distances. Hence, even for an average separation 
between heavy atoms, R, close to 2/ko, they approach each other to distances 
smaller than with a small tunneling probability P cx exp(— /3-\/ M/ m) <^ 1, 
where /3 ~ 1. It is then possible to extend Uj^^{R) (or U^j^{R)) to R < Kq^ 
in a way providing a proper molecule-molecule scattering phase shift in vacuum 
and verify that the phase diagram for the many-body system is not sensitive to 
the choice of this extension [68j . 

In Fig. 13 we display the zero-temperature phase diagram obtained by the 
Diffusion Monte Carlo method [SS]- Simulations were performed with 30 parti- 
cles and showed that the solid phase is a 2D triangular lattice. For the largest 
density it has been verified that using more particles has little effect on the 
results. 

For both 2x3 and 2x2 cases the (Lindemann) ratio 7 of the rms displacement 
of molecules to R on the transition lines ranges from 0.23 to 0.27. At low 
densities n the de Broglie wavelength of molecules is A '--^ jR ^ k^^, and 
U^^{R) (or U^^{R)) can be approximated by a hard-disk potential with the 
diameter equal to the 2D scattering length. Then, using the DMC results for 
hard-disk bosons [77j . we obtain the transition lines shown by solid curves in 
Fig. 13. At larger n, we have A < Kq ^ and using the harmonic expansion of 
U{R) around equlibrium positions in the crystal, we calculate the Lindemann 
ratio and select 7 for the best fit to the Monte Carlo data points (dashed curves 
in Fig. 13). 

4.3 Molecular superlattice in an optical lattice 

The mass ratio above 100, required for the observation of the crystalline order 
can be achieved in an optical lattice with a small filling factor for heavy atoms. 
Their effective mass in the lattice, M*, can be made very large, and the discussed 
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Figure 13: Diffusion Monte Carlo gas-crystal transition lines for 3D (triangles) 
and 2D (circles) motion of light atoms. Solid curves show the low-density hard- 
disk limit, and dashed curves the results of the harmonic approach (see text). 

solid phase should appear as a superlattice. There is no interplay between the 
superlattice order and the shape of the underlying optical lattice, in contrast to 
the recently studied solid and supersolid phases in a triangular lattice with the 
filling factor of order one [781 ES HQ] . The superlattice discussed in our review 
remains compressible and supports two branches of phonons. 

The gaseous and solid phases of weakly bound molecules in an optical lattice 
are metastable. As well as in the gas of such molecules in free space, the main 
decay channels are the relaxation of molecules into deep bound states and the 
formation of trimer states by one light and two heavy atoms. The relaxation 
into deeply bound states turns out to be rather slow, with a relaxation time 
exceeding 10 s even at 2D densities ~ 10^ cm~^ [68 . 

The most interesting is the formation of the trimer states. In an optical lat- 
tice the trimers are eigenstates of the Hamiltonian Hq = — (?i^/2M*) 2 ^R.i + 
e+(i?i2). In a deep lattice it is possible to neglect all higher bands and regard 
Ri as discrete lattice coordinates and A as the lattice Laplacian. Then, the 
fermionic nature of the heavy atoms prohibits them to be on the same lattice 
site. For a very large mass ratio A'h/m the kinetic energy term in i/p can be 
neglected, and the lowest trimer state has energy etr ~ where L is the 
lattice period. It consists of a pair of heavy atoms localized at neighboring sites 
and a light atom in the gerade state. Higher trimer states are formed by heavy 
atoms localized in sites separated by distances R > L. This picture breaks 
down at large R, where the spacing between trimer levels is comparable with 
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the tunneling energy 'f? jM^l? and the heavy atoms are delocahzed. 

In the many-body molecular system the scale of energies in Eq. (|76p is much 
smaller than |eo|. Thus, the formation of trimers in molecule-molecule "colli- 
sions" is energetically allowed only if the trimer binding energy is etr < 2eo. 
Since the lowest trimer energy in the optical lattice is e+(L), the trimer forma- 
tion requires the condition e+(L) < 2eo, which is equivalent to Kq"^ > 1.6L in 
the 2x3 case and > 1.25L in the 2x2 case. This means that for a sufR- 
ciently small molecular size or large lattice period L the formation of trimers is 
forbidden. 

At a larger molecular size or smaller L the trimer formation is possible. 
The formation rate has been calculated in Ref. [55] by using the hybrid Born- 
Oppenheimer approach, and here we only present the results and give their 
qualitative explanation. In order to form a bound trimer state two molecules 
have to tunnel towards each other at distances R ^ Hq^- This can be viewed as 
tunneling of particles with mass M^, in the repulsive potential Ues{R)- There- 
fore, the probability of approaching at interheavy separations where the trimer 
formation occurs, acquires a small factor exp(— J-^/mT/w) with J ~ 1, and 
so does the formation rate. Thus, one can suppress the trimer formation by 
increasing the ratio M*/m. On the other hand, for M*/m < 100 these peculiar 
bound states can be formed on the time scale r < 1 s. 

Note that the trimer states in an optical lattice, at least the lowest ones, are 
much more long-lived than in the gas phase. An intrinsic relaxational decay is 
strongly suppressed as it requires the two heavy atoms of the trimer to approach 
each other and occupy the same lattice sites. The trimer state can also decay 
when one of the heavy atoms of the trimer is approached by its own light atom 
and another light atom at light-heavy separations ^ Re- However, this decay 
channel turns out to be rather slow, with a decay time exceeding tens of seconds 
even at 2D densities n ~ 10^ cm~^ [68] . 

5 Concluding remarks and prospects 

The most distinguishing feature of weakly bound bosonic molecules formed of 
fermionic atoms is their remarkable coUisional stability, despite they are in the 
highest ro-vibrational state. As we mentioned in the Introduction, the life- 
time of such molecules can be of the order of seconds or even tens of seconds 
at densities ^ 10^^ cm~^, depending on the value of the two-body scattering 
length. This allows for interesting manipulations with these molecules. One of 
the ideas is related to reaching extremely low temperatures in a gas of fermionic 
atoms at a < and achieving the superfluid BCS regime. This regime has not 
been obtained so far because of difficulties with evaporative cooling of fermionic 
atoms due to Pauli blocking of their elastic collisions. The route to BCS may 
be the following. In the first stage, one arranges a deep evaporative cooling of 
the molecular Bose-condensed gas to temperatures of the order of the chemi- 
cal potential. Then one converts the molecular BEG into fermionic atoms by 
adiabatically changing the scattering length to negative values. This provides 
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an additional cooling, and the obtained atomic Fermi gas will have extremely 
low temperatures < lO^^Tp, where Tp is the Fermi temperature. The gas can 
then enter the superfluid BCS regime [5T]. Moreover, at such temperatures 
elastic collisions are suppressed by a very strong Pauli blocking and the thermal 
cloud is in the coUisionless regime. This is promising for identifying the BCS- 
paired state through the observation of collective oscillations or free expansion 
[8l[83[Ti[16l[84j. 

It will also be interesting to transfer weakly bound molecules of fermionic 
atoms to their ground (or less excited) ro-vibrational state. For molecules of 
bosonic atoms this has been done using two-photon spectroscopy [SFj'SS',^ and 
by magnetically tuned mixing of neighbouring molecular levels, which enables 
otherwise forbidden radiofrequency transitions |88| . Long lifetimes of weakly 
bound molecules of fermionic atoms at densities ~ 10^'^ cm"'^ may ensure an 
efficient production of ground state molecules compared to the case of more 
short-lived molecules of bosonic atoms. One could then extensively study the 
physics of molecular Bose-Einstein condensation. Moreover, the ground-state 
heteronuclear molecules have a relatively large permanent dipole moment and 
can be polarized by an electric field. This may be used to create a gas of 
dipoles interacting via anisotropic long-range forces, which drastically changes 
the physics of Bose-Einstein condensation (see, e.g., [55J and references therein). 

In the last few years, the observation of the Efimov effect was one of the 
important goals in the cold atom studies. As we discussed in Section 3, the 
Efimov trimers in the gas phase are short-lived and rather represent narrow res- 
onances. The Efimov effect then manifests itself in the log-periodic dependence 
of collision rates on the two-body scattering length. In particular, this is the 
case for the rate of three-body recombination of atoms [7^ and for the rate of 
trimer formation in molecule-molecule collisions [57]. In this sense, the trimer 
formation in gases of bosonic molecules consisting of heavy and light fermions 
(such as, for example, LiYb) attracts a great interest as the observation of the 
Efimov oscillations requires a much smaller change of the two-body scattering 
length (by a factor of 7 or 5) than in the case of identical bosons. 

Of particular interest are the trimer states of two heavy and one light fermion 
in an optical lattice. For 2D densities ~ 10* cm^^ the rate of the trimer forma- 
tion can be of the order of seconds, and these states can be detected optically. 
As we already mentioned in Section 4, the lattice trimers are long-lived, with 
a lifetime that can be of the order of tens of seconds. Thus, it is interesting to 
study to which extent these non-conventional states, in which the heavy atoms 
are localized in different sites and the light atom is delocalized in between them, 
can exhibit the Efimov effect. 

The creation of a superlattice of molecules in an optical lattice also looks 
feasible. A promising candidate is the ^Li-''*'K mixture as the Li atom may 
tunnel freely in a lattice while localizing the heavy K atoms to reach high mass 
ratios. A lattice with period 250 nm and K effective mass M*= 20M provide 
a tunneling rate ~ 10'^ s~^ sufficiently fast to let the crystal form. Near a 
Feshbach resonance, a value a — 500 nm gives a binding energy 300 nK, and 
lower temperatures should be reached in the gas. The parameters tik^'^ of 
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Fig. 13 arc then obtained at 2D densities in the range lO'^ — 10^ cm ^ easily 
reachable in experiments. 
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